Why stripes? Spontaneous formation of inhomogeneous structures due to elastic 

interactions 
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can argue that in systems with only short-range interac- 
tion, such as the conventional Hubbard or t—J models, 
one would rather expect a total phase separation into 
two phases: an antiferromagnetic insulator without any 
holes (n = 1) and another phase — metallic and probably 
ferromagnetic — containing all the holes; this state allows 
to gain the full kinetic energy of holes, and simultane- 
ously it costs a minimum of surface energy. 

To stabilize the stripe phase, one often invokes long- 
range Coulomb interaction which prevents the large-scale 
phase separation gjlfj. However, for the strong long- 
range repulsion we would rather expect not stripes but a 
sort of Wigner crystal, in which the doped carriers keep as 
far away from each other as possible. In this case, stripes, 
and even more so two-dimensional sheets observed e.g. 
in manganites |5|,||, are definitely less favourable. One 
then argues that stripes would be stable in the interme- 
diate case. And although this scenario is of course not 
excluded, it seems that some important physical factors 
are missing from such a picture. ^ 

We argue that there is indeed one such factor, which 
exists in all real systems and which can, in principle, 
stabilize the stripe and sheet phases. This factor is the 
long-range elastic interaction. We consider below the 
insulating materials (nickelates, overdoped manganites) 
with the so called filled stripes (one hole per site); the 
situation in cuprates with presumably half-filled stripes 
may be more complicated, although we believe that the 
factors invoked below can also play an important role 
there. 

When we dope a stoichiometric system, e.g. substitute 
Ca for La in LaMnC>3 , we create Mn 4+ ions in the Mn 3+ 
matrix. Or, in the overdoped Lai_ x Ca x Mn03 with x > 
0.5 (this is the situation, in which the stripe, or rather 
sheet phases were observed in manganites ) , one can 
speak of a certain amount of Mn 3+ ions in the Mn 4+ 
matrix. 

Classically, Mn 3+ ions differ from Mn 4+ (or Ni 2+ from 
Ni 3+ in nickelates |Rpl) not only by their charge, but 
also by their ionic radius. This factor is usually ignored 
in the conventional treatment of stripes. In contrast, we 



We argue that elastic interactions between ions in differ- 
ent valence states can play an essential role in stabilization of 
stripes (or 2D "sheets") in doped oxides. These interactions 
are in general long-range and anisotropic (attractive in cer- 
tain directions and repulsive in others). This can naturally 
give rise to stripe-like structures in insulating materials. We 
illustrate this general idea with certain specific examples and 
show that the situation can be described by the Ising model 
with anisotropic interactions. The case of anisotropic impu- 
rities, relevant e.g. for manganites, is also briefly discussed. 
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The notion of stripes became recently one of the most 
important concepts in the physics of doped Mott insula- 
tors. There are many theoretical and experimental 
[^|-^| indications that they indeed arise in doped transi- 
tion metal oxides; there exist also serious arguments that 
they may play an important role in high-temperature su- 
perconductivity [p]j9|JTo[| . 

Despite very intense studies in this field, the origin of 
stripes in real materials is far from being well understood. 
Following the pioneering work JlJ the main attention is 
paid to the purely electronic and magnetic mechanisms of 
stripe formation However the existence of stripes 

in this approach is still controversial . 

When discussing the origin of stripes in doped oxides, 
one usually starts with Mott insulators with an integer 
number of electrons per site forming local magnetic mo- 
ments. Generically their ground state is antiferromag- 
netic. In this case, one argues that with doping a su- 
perstructure may develop consisting of antiferromagnetic 
domains, the doped charge carriers being localised at the 
domain walls. This reflects the general tendency: insta- 
bility of a homogeneous state in doped strongly corre- 
lated electron systems toward phase separation [[IT|-[h4[ . 
Stripe phases are possible manifestations of this ten- 
dency. 

Why are stripe phases better than any other forms of 
phase separation? This is the main question, which, in 
our opinion, still did not find a satisfactory answer. One 
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2 Note also that in nickelates and manganites, in contrast to 
cuprates, the metal-centered stripes would not correspond to 
the antiferromagnetic domain walls inherent to the conven- 
tional theories fllr 



1 



pay the main attention here just to this factor (which, 
of course, should be considered alongside with the other 
factors, usually accounted for JT| |^,|o|,p~0|] ) . 

As is known from the theory of elasticity, when we 
create a lattice distortion (in simple case — a dilatation, 
e.g. by cutting out a small sphere of radius A and sub- 
stituting it by the sphere of radius A' / A), we also 
create a field of lattice strains, which is long-range and, 
in general, anisotropic |l5|,[l6| . This leads to a long-range 
interaction between impurities (~ 1/i? 3 ), which depends 
on the elastic constants of the medium and on the shape 
of impurities. 

The most important point is that typically this inter- 
action resembles a quadrupole-quadrupole interaction, 
being repulsive in certain directions but attractive in 
others. This long-range attraction opens a possibility 
to form regular structures — clusters of impurities form- 
ing one-dimensional objects (needles or stripes) or two- 
dimensional structures (lamellae or sheets); this can be 
one of the mechanisms stabilizing stripe or sheet phases. 

The interaction of dilatation impurities in crystals van- 
ishes in isotropic media, but it is nonzero in crystals. 
Thus, in crystals with cubic anisotropy it is given by |l5| ] 

T(n) 

V{^) = -CdQ x Q 2 r^^ (1) 

|?' — r \ A 

where C is a constant of the order of unity, f — r 1 = 
\f — r*| • ft, Qi and Qi are the "strengths" of impurities 
(Qi ~ ( v i ~ v o) where Vi is the volume of the impurity 
and Vo is the corresponding volume of the matrix), and 

d = en - C12 - 2c44, (2) 

where Cy are the elastic moduli of the crystal. The angu- 
lar dependence of the interaction ([!]) is determined by a 
function of the direction cosines of the vector R = ?\ — f^'- 

T{ri)=ni+n 4 y +ni-\. (3) 

From (0)-(§) we can see that the interaction between im- 
purities is long-range (~ 1/i? 3 ) and has different signs in 
different directions. Thus, for the situation with d > 
it is attractive along [100], [010] and [001] directions 
(r([100]) = |) and is repulsive along [110], [011] etc. 
(r([110]) = —jo) and along the cube diagonals [111]: 
r([lll]) = —jg, and vice versa for d < 0. Thus, if we put 
a few dilatation impurities in a cubic crystal, they will at- 
tract each other along certain directions (e.g. [100], [010], 
[010]) and repel along others, thus causing a formation 
of inhomogeneous structures (stripes, sheets etc.)P| 



3 Note also that in contrast to the Coulomb forces, these 
elastic interactions are not screened (except by the "mirror 
forces" due to the surface of the crystal Pa]). 



In real systems (not the weakly anisotropic crystals 
considered above as an example) , the ratios of the inter- 
action constants between different neighbors may differ 
from those given above; thus, e.g., in bec iron the in- 
teraction with [011] neighbor is stronger than with [001] 
one But again, as in the example treated above, 

the interactions have different signs in different direc- 
tions. Thus, we can consider the situation of a cubic 
lattice with the nearest-neighbor (nn) and next-nearest- 
neighbor (nnn) interaction of opposite signs as a typical 
one, treating the ratio of these couplings as an arbitrary 
parameter. This naturally leads to the lattice gas, or the 
Ising model of the type 

H = — J <Ji(Tj + J' ai<Tj + const (4) 

(ij)-nn (ij)-nnn 

where (jj = ±1; the density at site i is rii = \{l+(Xi), and 
the model has to be considered for fixed density n = (n*) 
(or fixed "magnetisation"). We treat below the case of a 
cubic lattice. 

For J, J' > (which corresponds to d > in Eq. ([j])) 
we have the following solutions: 

1. Totally phase separated state, in which all the "im- 
purities" (or electrons) form one big cluster. The energy 
of this state (per impurity) in a simple cubic lattice is 
£ph.sc P . = -Jz nn /2 + J'z nnn /2 = -3 J + 6 J'. Here z are 
the corresponding numbers of nearest and next nearest 
neighbours. 

2. Occupied sites forming 2D sheets parallel to (xy), 
(xz) or (yz) planes, such that the distance between them 
exceeds one lattice spacing. The energy of such state is 
-Eshcct = ~2 J + 2 J'. 

3. The II? stripes in x, y or z directions, which do 
not cross or approach one another to one lattice spacing. 
The energy of such a state is -E s tripe = — J ■ 

One can easily show that all the other possible struc- 
tures have higher energy. 

By comparing the energies of these states, we see 
that the totally phase separated state will be stable if 
J 1 1 J < \; for i < J'/J < \ the 2D "sheets" are stable, 
and for J'/J > i the ID stripes would give a minimum of 
the energy. Thus, we see that the sheet or stripe phases 
appear quite naturally if the dominant interaction is the 
elastic interaction between impurities which is intrinsi- 
cally anisotropic and long-range. For d < 0, we similarly 
obtain, instead of the vertical, the diagonal stripes, but 
in this case, we should also include the interaction with 
third neighbors. 

For low electron (or impurity) density in 3D case the 
lamellae phase is formed by parallel sheets; the distance 
between them in model (|4|) with only nn and nnn in- 
teractions is arbitrary, provided only that it exceeds the 
lattice spacing. Longer-range interactions will make this 
arrangement regular, but for weak interaction of this kind 
the discommensurations could be easily formed. 

The situation is a bit more complicated for the stripe 
phase. In layered materials, it will be similar to that 
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with the sheets in 3-D crystals (parallel equally spaced 
stripes). In 3D crystals there may in principle exist mu- 
tually perpendicular stripes. Again, one may expect that 
longer-range interaction can make these stripes parallel; 
this question, however, goes beyond the scope of our 
treatment. Thus, we see that the elastic interaction of 
the "spherical" impurities in crystals may quite naturally 
lead to the formation of inhomogeneous structures — ID 
stripes, 2D sheets or 3D phase-separated clusters, which 
can be described by the Ising-like model (||) [|l7|. Note 
also that this process does not require the atomic diffu- 
sion: it is realized by the electron hopping between, say, 
Mn 3+ and Mn i+ (with the corresponding lattice relax- 
ation). 

The problem of formation of regular structures, when 
we have two components with different atomic volumes, 
is well-known in the physics of segregating alloys. Theo- 
retical studies (see e.g. |l6)) have shown that the shape 
of inclusions of one phase in another is to a large extent 
determined by the elastic strains: the shape is such as to 
minimize the total elastic energy. As is shown in the cor- 
responding studies, the best way is to create the second 
phase in the form of infinitely thin layers having certain 
specific orientation in given matrix. In the correspond- 
ing alloys it gives rise to the so called Guinier-Preston 
zones flgjl^] : inclusions of the second phase appear as a 
regular array of parallel platelets. One can speculate that 
the physics of the stabilization of stripe phases (actually 
sheets) in Lai-^Ca^MnOs at x = 0.67, 0.75, etc. is 
the same, and that it is just these long-range anisotropic 
elastic forces which stabilize such structures. This is sup- 
ported by the fact that this stripe charge ordering is ob- 
tained in the insulating materials at temperatures above 
the eventual magnetic ordering; thus, the often invoked 
magnetic mechanisms play here a minor role, if any. The 
inclusion of the electron hopping will, of course, tend to 
destabilize these ordered structures, which could survive 
if the hopping integrals are not too large. 

The situation may be different (and actually much 
richer) if the "impurities" are not simple dilatation cen- 
ters ("spheres"), but are anisotropic ("ellipsoids"). This 
is the typical situation with Jahn- Teller ions, e.g. Mn 3+ 
"impurities" in overdoped manganites Ri-^Ca^MnOa (R 
= La, Pr, Nd) for x > 0.5. The interaction between such 
impurities also decays as I /R 3 and has different signs de- 
pending on both the relative position of such centres in a 
crystal and on the orientation of corresponding orbitals 
(i.e. local distortions) (2^^l[]. Thus, for example, for two 
quadrupolar ions along the z-axis with the electron den- 
sity elongated parallel to x, y, or z axes, as is the case of 
e 9 -ions_Jike Mn 3+ one obtains from the general expres- 
sions 



the interaction in the form 
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(en + c 44 ) 
8tt(cii + 2 Cii )R 3 
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where en and C44 are elastic moduli, and a aa is a stress 
tensor such that e.g. for the center with the occupied 



orbital 3z 2 — r 2 we have a zz = 1, <r. 



-1/2 (and 



corresponding expressions for z — > x, y). 

From this expression one can see that e.g. two 3z 2 — r 2 
orbitals along z-axis strongly repel one another, whereas 



3z l 



r2 and 3x 2 



(or 3y 



') along the same di- 



rection attract. This can quite naturally explain the 
well-known orbital structure of undoped LaMnOs (cf. 
p2j), and also a special stability of the Mn 3+ -planes, 
with the corresponding orbital ordering, in lightly doped 
manganites, e.g. Lai-^Sr^MnOs, x ~ 1/8 [^5| and in 
Pri.zCa^MnOa, x ~ 1/4 j26|. With the assumption 
of a checkerboard charge ordering it also gives the or- 
bital ordering of CE type observed in R 1 _ 2 ,Ca :! ,Mn03 at 
x = 0.5 ]23| , p4] |. Indeed, in both cases the orbitals of the 
neighbors along x and y directions are mutually orthog- 
onal (e.g. 3x 2 — r 2 and 3y 2 — r 2 orbitals); according to 
(^|) they have an attractive interaction which stabilizes 
these structures. One can also show that the interaction 
along diagonals ([110] direction) is attractive for the par- 

but repulsive for the 
r 2 at 



allel orbitals 3x — r or 3y 
perpendicular ones (3x 2 



at one site and 3y 
the other); this gives extra stability to the observed or- 
bital ordering in undoped LaMnOs, and these diagonal 
interactions nearly cancel for the CE-type phase in the 
x = 0.5 system. The same factors may favor single stripes 
(" Wigner crystal" ) || as compared to the bi-stripes (H) . 

Summarizing, we suggest that the elastic interactions, 
always present in systems with mixed valence, e.g. in 
doped Mott insulators, may play an important role in 
stabilizing particular types of inhomogeneous structures, 
such as stripes or two-dimensional sheets. These inter- 
actions decay rather slowly (~ R~ 3 ) and are anisotropic, 
being repulsive in some directions but attractive in oth- 
ers. Thus, the interaction of dilatation impurities in 
weakly anisotropic cubic crystals has different signs along 
[100], [010], and [001] directions and along face or body 
diagonals. Such an interaction may quite naturally lead 
to the formation of stripes or sheets in manganites or 
nickelates. Similar interactions in the case of anisotropic 
(Jahn- Teller) ions can produce analogous inhomogeneous 
structures with specific orbital ordering. The importance 
of lattice effect in the stripe formation is supported by 
the observation of a large isotope effect on the stripe or- 
dering temperature [ p7|p8| ]. Thus, we can conclude that 
the elastic or, more generally, electron-lattice interac- 
tions may be very important in providing a mechanism 
for stripe formation in doped systems. 
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